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Abstract 

Unstable homoepitaxy on rough substrates is treated within a hnear contin- 
uum theory. The time dependence of the surface width W{t) is governed by 
three length scales: The characteristic scale Iq of the substrate roughness, the 
terrace size lo and the Ehrlich-Schwoebel length Ies- Ies Id (weak step 



edge barriers) and /q ^ ~ IdV^d/Tes, then W(t) displays a minimum at 

a coverage 6'inin ~ {Id/Ies)'^, where the initial surface width is reduced by 

a factor lo/lm- The role of deposition and diffusion noise is analyzed. The 

results are applied to recent experiments on the growth of InAs buffer layers 

[M.F. Gyure et ai, Phys. Rev. Lett. 81, 4931 (1998)]. The overah features 

of the observed roughness evolution are captured by the linear theory, but the 

detailed time dependence shows distinct deviations which suggest a significant 

influence of nonlinearities. 
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I. INTRODUCTION 



A high symmetry crystal surface growing epitaxially from a molecular beam can become 
unstable towards the formation of pyramidal mounds if the mass transport between different 
atomic layers is reduced by additional energy barriers at step edgesllll. Over the last few 
years, this phenomenon has been observed for a wide range of metal and semiconductor 
surfaces, and a considerable body of theoretical work has been devoted to the description 
of the asymptotic (late time) evolution of the surface morphologyiHl. In the early time 
regime continuum theory predicts an exponential growth of the surface modulations. For 
this reason the precise initial state of the surface has commonly been disregarded, since the 
exponential instability should rapidly wash out the details of the substrate roughness. 

In a recent paper! Gyure, Zinck, Ratsch and Vvedensky (GZRV) presented experimental 
and numerical results for the early time development of unstable homoepitaxy from a rough 
substrate, which show a more complex scenario: It was observed that the competition 
between smoothening of the initial roughness and the instability associated with the incipient 
mound structure can lead to a minimum in the total surface width. A similar effect was 
predicted previously in the context of noise-induced roughening^, and related experimental 
observations have been reported both for thin metal filmsi and semiconductor multilayersi. 
Qualitatively, the minimum originates from the wavelength dependence of smoothing and 
(deterministic or stochastic) roughening rates: If the roughness spectrum of the substrate 
has sufficient weight at short wavelengths, which are efficiently smoothened by capillarity 
effect Jl^, then the decrease of the substrate contribution to the surface width can temporarily 
dominate the long wavelength roughening induced by growth. 

The possibility to minimize the surface roughness by an appropriate choice of the buffer 
layer thickness and other growth parameters is of obvious interest in applications. In this 
paper we develop a quantitative theory of unstable growth on rough substrates, which allows 
us to determine the conditions under which a minimum occurs, and to estimate the layer 
thickness of minimal roughness in terms of microscopic length scales and parameters, such 
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as the in-layer and interlayer diffusion barriers. Our starting point is the observation of 
GZRV that the time evolution of the roughness spectrum appears to be well described by 
the linearized continuum evolution equation for the surface. By incorporating various kinds 
of nois into the linear theory we can compare the influence of stochastic and deterministic 
roughening, and obtain a unified description of both cases. A critical discussion of our results 
in relation to the experiments of GZRV will be presented at the end of the paper. 



II. LINEARIZED CONTINUUM THEORY 

The standard phenomenological evolution equation for the continuous surface profile 
H{r, t) is of the formEH 

dtH + V ■:i = F (1) 

where the surface current J incorporates both a growth-induced destabilizing 
contributionHEl and a stabilizing term originating in capillarity0, and F denotes the de- 
position flux, which will be assumed constant for the time being. Small fluctuations hijc, t) 
around the fiat singular surface H = Ft then satisfy the linear equation 

d^h = -aV^h - K{V^fh (2) 

with positive coefficients a, k representing deposition (a) and smoothening (k), respectively, 
whose relation to the growth parameters will be explained below. 

The substrate roughness is incorporated through a spatial roughness spectrum 
(|/i(k, 0)p) = S{k,Q) = So{k), where h{]i,t) is the Fourier transform of h{T,t) and k = |k|. 
Under the linear equation (0) the roughness spectrum evolves as 

S{k, t) = So{k) exp[2{ak^ - Kk^)t], (3) 

which implies that fluctuations with wavenumbers k > kc = ^Ja/ k are damped while those 
with k < kc are exponentially amplified. The surface width W{t) is obtained by summing 
over all wavenumbers, 
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W\t) = 211 r dk A;^o( A;) e^^"'^' (4) 
Jo 

Motivated by the experimental data shown in Figure 4 of GZRV we choose a white noise 
roughness spectrum, 

■ (5) 

: else, 

where the small scale cutoff Iq is required for a finite value Wq = W{0) of the initial surface 
width. Taking the time derivative of @ and evaluating it at t = 0, we find that the surface 
width shows an initial decrease if {Im/loY > 12, where = 2^^^J2^^/a is the wavelengthil 
of those fluctuations which are maximally amplified by the linear equation (0). Thus the 
condition for a nonmonotonic time dependence of the surface width is that the length scale 
characterizing the substrate roughness, Iq, is much smaller than the typical scale Im of the 
emerging mounds. This result holds also for more general initial roughness spectra, e.g. 
S'o(fc) = Ak~^ with p < 2 and a small scale cutoff Iq. For substrates whose roughness is 
dominated by long wavelength fluctuations, in the sense that 5*0 ~ k'^ with p > 2, a large 
scale cutoff is needed and the time derivative dW/dt\t=o does usually not exist. 
In the following we take Iq <^ Im- Then Eq.(^) reduces to the scaling form 

W\t) = W'Mlrr;)H{t/T), (6) 

where 1/r = a^/4K is the amplification rate of the maximally unstable fluctuations and the 
scaling function is 

<l>(x) = e'^yVa^ll + erf(v^)] (7) 

with erf (s) = /q* exp(— The width attains its minimum at a time tmin ~ 0.18 r, 

where it has been reduced by a factor 

W{t^^)/W, ^ 3.7 [l^/U. (8) 

Since the factor 1 + erf (-\/2x) in (|^ only varies between 1 and 2, the scaling function 

is essentially the product of a decaying power and an exponentially increasing factor. The 
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power law for small x reflects the particular smoothening mechanism (capillarity-driven 
surface diffusion) and its general formS is given by ([T8|) below. For finite lo/lm the power 
law sets in for times t > to with to ~ (^o/^m)^T- 

To relate the behavior of W{t) to microscopic parameters we need to express the coeffi- 
cients a and k of (|^) in terms of the two length scales governing unstable homoepitaxyii: 
The typical terrace size^l and the Ehrlich-Schwoebel-lengthlii 

Ies = aii{D/D' - 1) = a,|(e^^/'=^^ - 1) (9) 

defined in terms of the in-layer lattice spacing ay, the in-layer (interlayer) surface diffusion 
constant D [D') and the step edge barrier /\E. Comparison of the two length scales allows 
to distinguish conditions of strong {Ies ^ Id) and weak {Ies ^ Id) step edge barriers; in 
the first case a ~ Ftj-^, in the second a ~ FI^Ies- The coefficient k, is traditionally associ- 
ated with near-equilibrium surface diffusionS, however under far-from-equilibrium growth 
conditions the dominant contribution to k, is believed to arise from the random nucleation 
process^!. The expression k Fljj is then suggested by dimensional analysis^ and scaling 
arguments^. It leads to a consistent picturei in the sense that I^k, and r in the 

strong barrier case, which implies that mounds develop on the submonolayer islands already 
during the growth of the first few layers ("wedding cake" regime@'0). In the weak barrier 
case we find 



Ira ~ Id\JId/Ies and r ~ F-^h/hs)^. (10) 
The minimum in the surface width thus occurs at a coverage 

^min ~ (Id/Ies)^ > 1 (11) 

which corresponds, not surprisingly, to the coverage where mounds first become visible for 
growth from a smooth substratei'0. Similarly the coverage 6q = Fto at which the scaling 
form (^ for the width begins to hold is of the order of Oq ~ {Iq/Id)^ independent of Ies 
(provided Ies ^ Id)- 



To apply these considerations to the experiment on InAs growth of GZRV, we first 
need to check the condition Iq -C Im- /,From Figure 4 of the paperl we estimate that 
Wo/W{t^in) ~ 4. Comparing this to the theoretical prediction (|) we find Iq ~ 0.07 x 1^ 
and Im/^o ^ 1 is true. This is in contrast to the kinetic Monte Carlo simulations of GZRV, 
where W()/W{t^i^) ~ 1.1. The instability length in the experiment is ~ 1.0 /im, which 
yields /q ~ 70 nm for the small scale cutoff of the substrate roughness. This is consistent 
with the initial roughness spectrum in Figure 4 of GZRV, which is constant at least down 
to a length scale of 300 nm. The minimum width is attained at a film thickness of about 
0.57 /im. Using a\\ ~ 6A and a bilayer thickness a± ^ 3A, we therefore estimate that 
^min ~ 1900 and ~ 1700, and hence /_BS'/a|| ~ 6 and /d/«|| ~ 250. At the experimental 

temperature of 500° C this implies a step edge barrier AE of the order of 0.1 eV, comparable 
to estimatesi'0 for GaAs. 



III. NOISE EFFECTS 

Next we include a noise term ri{r,t) in Eq. (j^). The different sources of noise, the 
individual events of deposition ("shot noise") and diffusion, enter the noise correlator with 
different dependenceE'lil on the wavenumber k. We write it in the form 

R{k) = (r/(k, t)r/(-k, t)) = Rs + RBk' (12) 

with Rs and Ro denoting the strength of deposition and diffusion noise, respectively. In the 
linear model with noise the roughness spectrum S{k,t) then contains a part reflecting the 
history of the noise, as well as the deterministic evolution of the initial roughness treated 
above. The full expression reads 

S{k,t) = Sdetik,t) + Snoiscik,t) = (13) 



So{k) e2("^'^ + 



2{ak^ - Kk^) 

Unlike the deterministic mechanism, the noise increases the amplitude of the spectrum for 
every wavelength, i.e. dtSnoisc{k,t) > for all k. We shall now examine whether under the 



experimental conditions of GZRV noise substantially contributes to the surface width. 

The deposited particle flux can be seen as a Poisson process with intensity F, so Rs = 
a±a'jF. Shot noise thus contributes to the total width by 

Wsity = FTa^{a\\/l^f^>{t/T), (14) 

with \E''(x) = for the choice (|^) of 5*0 (/c). Using (0) we see that (p^if) can be can be 
ignored against Eq. if 

(W^oK)'(/o/«||)' » {Id/Ies?. (15) 

With our estimate /q ~ 70 nm this condition is satisfied in the experiment. A different 
interpretation of Eq. (|l^) will be given below. 

The diffusion noise strength is given by the average rate of adatom jumps on the surfaced, 

so 

Ri^^piD^llF (16) 

where we have used the estimate pi ~ Fl'j^/D for the adatom densitylii. Diffusion noise 
thus becomes more important than shot noise for k>7r/l£), whereas it can be neglected for 
long wavelengths. For large k we can approximate the contribution of diffusion noise in (|T3|) 
by Rj^k'^/ink'^) which enters the total width as WDit)"^ = /|)//^-Fr log(/£)/a||). This means 
that roughly Wij{t)^ ^ (^D/^m)^Ws(^)^ ^ W^s(t)^ because > in the weak barrier 
regime. In particular, at the time when the width minimum is attained, diffusion noise can 
be neglected against shot noise and for the experiment of GZRV Eq. remains valid. 

It was mentioned already that due to noise a minimum in the surface width may occur 
even in the absence of step edge barrier^l. For completeness we provide here a simple analysis 
for the most general linear Langevin equation of kinetic roughening, 

d^h = -{-uVy/^h + r] (17) 

where z = 2 and z = 4 correspond to evaporation-condensation and surface diffusion domi- 
nated relaxation, respectively^, i/ > is a constant and and r] is the deposition noise. Odd 



or noninteger values of z describe nonlocal relaxation mechanisms and can be treated on the 
same footingl. The linearity of ( p!?]) implies that the substrate contribution and the growth 
induced contribution to the roughness can be separate d@. The substrate contribution is 
found to decay according to 

vr,,b(t) ^ w^ihimYi'' (18) 

for a (i-dimensional surface, for times such that the correlation length of the growth-induced 
roughness ^(t) exceeds /q (otherwise PVsub ~ Wq)- In terms of physical quantities the corre- 
lation length can be written a s0 i{t) ^ Id^^'"'- To determine the coverage Qrmn of minimal 
surface width for purely stochastic roughening, the substrate contribution (|18D should be 
compared to the growth induced roughness0 

W^growth^ax(0/^~)(^-'^)/2^ (19) 

where Q is the coverage at which the width becomes of the order of a_L and thus lattice 
effects (such as temporal oscillations of the step density) die out; the expression ([T9|) holds 
for ^ > ^. With the estimateBi ~d ~ {l^j a\\)^'^!^^~'^^ we obtain 

^min~ (W«±)'(^o/a||)' (20) 

independent of z. Comparing Eqs. (|20| ) and (|Tl|) we have thus recovered the crossover 
condition (|T5|) between deterministic instability and stochastic roughening from the opposite 
side. To neglect the growth instability in Eq.([T7|) is no longer justified when the minimum 
width coverage 6'min predicted by the deterministic theory (Eq.(|Tl|)) is smaller than 6'mm in 
Eq.(©. 

IV. DISCUSSION AND CONCLUSION 

The main results of this paper are Eqs. (^ and (0), which express the time dependent 
surface roughness in terms of the characteristic length and time scales of the problem - the 
substrate roughness scale Zq; the incipient mound size Im-, and the linear growth time r -, 
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the latter two of which are, in turn, related to the microscopic growth parameters through 
Eg. (p!OD ■ For the experiment of GZRV, the measured values of 1^ and r were seen to imply 
reasonable numbers for the microscopic lengths Id and Ies, and for the step edge barrier 
AE. 

It is then natural to ask to what extent the linear theory can be used to quantitatively de- 
scribe the experimentally observed roughness evolution, beyond providing consistent order- 
of-magnitude estimates. Inspection of the data for the surface roughness depicted in the 
inset of Figure 4 of GZRV quickly leads to the conclusion that, despite a similar overall 
appearance, the shape of W{t) is not well reproduced by our scaling functions (|) and (|14|). 
In fact the data for the structure factor in Figure 4 show a qualitative feature which the 
linear theory is unable to explain: It is an immediate consequence of Eq.([13|) that S{k,t) 
is a monotonic function of t (increasing or decreasing) for any k; in contrast, the measured 
structure factor shows a nonmonotonic dependence on film thickness for k > kc- 

This prompts the question whether the use of the linearized theory is really justified 
under the experimental conditions. Nonlinear terms in the surface current J in Eq.(|I]) are 
expectedi00 to matter when the surface slope \Vh\ becomes comparable to a±/l£). Since 
typical slope values of the initial surface profile are of the order of Wq/Iq, the condition for 
the validity of the linear theory is 

Wo/a^ < lo/h. (21) 

In the experiment of GZRV Wo/a±_ ^ 3 and, from the estimates presented above, Iq/Id ~ 
0.5; thus the condition (^) is (weakly) violated. The analogy to phase ordering kinetic^l 
suggests that the early time evolution may be qualitatively altered when nonlinearities are 
important!^. This seems like an interesting problem for future study. 
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